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Transgenic Bacillus thuringiensis (Bt) crops play an increasing role in pest control, and resistance 
management is a major issue in large-scale cultivation of Bt crops. The fitness cost of resistance in targeted 
pests is considered to be one of the main factors delaying resistance when using the refuge strategy. By 
comparing 10 resistant Helicoverpa armigera (Hubner) strains, showing various resistance levels to Bt toxin 
(Cryl Ac), to a susceptible strain, we showed an increasing fitness cost corresponding with increasing levels 
of resistance. The relationship between overall fitness cost C and the resistance ratio R r could be described by 
C = 24.47/(1 + exp([1.57 - Log 10 J? r ]/0.2)). This model predicted that the maximum overall fitness cost 
would be —24% (±5.22) in the strains with the highest resistance level. The overall fitness cost was closely 
linked to egg hatching rate, fecundity, emergence rate, larval survival rate, and developmental duration of 
adults. Among fitness components measured, fecundity was the most sensitive trait linked to the resistance 
selection. To integrate the results into simulation models would be valuable in evaluating how variation in 
fitness cost may influence the development of resistance in pest populations, thus helping to develop 
enhanced refuge strategies. 



Transgenic plants producing the Bacillus thuringiensis toxins (Bt), which primarily target lepidopteran pests, 
have been planted over widespread areas for almost two decades. Bt corn and Bt cotton, for instance, covered 
69 million hectares globally in 20 12 1,2 . Bt crops have led to a drastic drop in insecticide use as well as steadily 
increasing biocontrol services provided by generalist predators 3 5 . Widespread adoption of Bt crops has also 
benefited neighboring non-Bt crops through suppression of major insect pest populations and conservation of 
natural enemies 3 5 . Nevertheless, sustainability of Bt crops may be compromised due to the development of 
resistance by pests, e.g. field-evolved resistance has occurred in populations of several major targeted pests 7-8 
years after Bt crops were first implemented 6,7 . Various theoretical mathematical models have suggested that the 
refuge strategy 8 , if widely adopted, would be effective in delaying insect resistance to Bt crops 6,7 . Analyses of more 
than a decade of global resistant monitoring data matched with these predictions as most targeted pests remained 
susceptible to Bt crops when the refuge strategy was used 6,7 . 

The basic theory underlying the refuge strategy is to reduce resistance heritability by providing numerous 
susceptible moths in refuges to mate with the resistant moths from Bt crops 6-9 . To delay pest resistance, the main 
method is to promote survival of susceptible pests by planting refuges of non-Bt host plants near Bt crops 6 8 . 
Fitness costs associated with Bt resistance, i.e. lower fitness in resistant insects than in susceptible ones in the 
absence of Bt toxins, can further delay the development of resistance by selecting against Bt-resistant genotypes in 
refuges, and this strategy can even reverse resistance development 10 . 

Resistance to Bt is controlled by multiple genes, and multiple resistance alleles with different types of mutations 
have been found on the same locus 10,11 . If each mutation may cause a cost, and the costs may accumulate with 
more genes or alleles involved in conferring a higher resistance level. Thus, fitness costs are expected to increase 
steadily with the development of increased resistance. Although fitness costs have been identified in various 
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Figure 1 | Flowchart of selection of 1 1 strains of H. armigera from three field collected populations in the laboratory. The names of the strains 
were shown in rectangle frames, numbers of generations were shown in the diamond frames, doses of Bt toxins (mg CrylAc/L artificial diet) imposed on 
H. armigera larvae were shown on the arrows. The last generation of each strains were tested for fitness cost in the present study. 



targeted pest populations, with the higher costs found linked to 
higher resistance levels to Bt 10,12 , successive changes in fitness cost 
due to development of resistance to Bt have not been fully 
demonstrated. 

Fitness costs associated with Bt resistance are often detected by (i) 
comparing one fitness component, e.g. survival, development rate or 
weight, between Bt-resistant and Bt-susceptible strains in the absence 
of Bt toxins, or by (ii) evaluating the stability of resistance to Bt 
toxin(s) in heterogeneous laboratory strains reared on non-Bt 
plants 10 . The former method provides insight into which fitness 
components are affected by resistance, and the latter method mea- 
sures the cumulative effect of all fitness components over many 
generations 10 . However, these two methods do not enable clear iden- 
tification of actual fitness component(s) that may respond to selec- 
tion events for resistance to Bt toxins. Moreover, it leaves unknown 
how to compare the various costs detected with the two different 
methods. 

In the present study, ten resistant strains of H. armigera represent- 
ing different stages of resistance development to Cryl Ac were used to 
compare fitness cost with a susceptible strain to (1) depict how the 
fitness cost changes through resistance development and (2) identify 
sensitive fitness components related to overall fitness cost changes 
and the increase of resistance level through resistance development 
in H. armiger to CrylAc. 

Results 

Eleven H. armigera strains were established by screening field- 
collected populations with different dose of CrylAc toxins under 
laboratory conditions (Fig. 1). When compared with the susceptible 
strain 96S, the other ten strains showed different resistance levels 
against CrylAc toxin with the BtR showing the highest resistance 
level (3536.5-fold) and the other strains showing intermediated res- 
istance levels (Table SI). The intrinsic rate of population increase, r m , 
estimated on non-Bt diet, was the highest in the susceptible strain 
96S, and dropped in the other resistant strains showing various over- 
all fitness costs (Table SI). This enabled the assessment of a possible 
relationship between overall fitness cost and level of resistance to 
CrylAc toxin in H. armigera. 



With linear and non-linear regression analysis, a three-parameter 
logistic model C = 24.47/(1 + exp([1.57 - Log lo K r ]/0.2)), with the 
lowest AIC (Akaike's Information Criterion) value of 81.77, was 
established to describe the relationship between overall fitness cost 
and resistance level (Fig. 2). This logistic model predicted that the 
overall fitness cost would increase to a limit of 24.47% (±5.22) in the 
highest CrylAc-resistant H. armigera strains. 

Thirteen fitness components were estimated for each strains with 
life tables. Comparing with the susceptible strain 96S, fitness com- 
ponent costs were often observed in resistant strains, i.e. they had 
lower survival rates, lower copulation rate, lower fecundity and 
longer developmental durations (Fig. 3 and 4), with exceptions for 
sex ratio and developmental duration of male and female adults, 
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Figure 2 | The increase of fitness cost with development of resistance in 
H. armigera to CrylAc. The dotted line indicates the linear tendency 
C = 0.09 + 8.71 Log 10 R r (p = 0.023, AIC = 82.50), while the solid line 
indicates the logistical tendency C = 24.47/(1 + exp( [ 1 .57 - Log 10 R r ] /0.2)) 
(AIC = 81.77). 
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Figure 3 | Survival rate from neonates to the 6 th instar larvae (A) survival from 6th instar larvae to pupa (B), emergence rate (C), copulation rate (D), 
hatching rate (E) and proportion of female adults (F) for each strains. The resistance level increased through the strains 96S to BtR. Bars indicate 
the 95% Credible Interval. Histograms bearing different letters are significantly different at P < 0.05. The statistical analysis is based on Bayesian methods 
and Monte Carlo approximations (see Table S2, S3-S8 for detailed statistical results). 



where no significant differences were found among the 1 1 strains on 
non-Bt diet (Fig. 3F and 4D & E). However, the variation of 10 fitness 
components showed various complicated trends with the increase of 
resistance to Cryl Ac (Fig. 3A-E and 4A-C & F-G, Fig. S1KA-KJ, Fig. 
S2KA-KJ & LA-LJ) or the overall fitness cost (Fig. S1LA-LJ) and thus 
were analyzed with multiple regressions to find the most influential 
fitness components. 

A multiple regression model C = — 1.591 + 0.03CW* + 0.258Cp ec 
+ 0.618Q;, + 0.618CR, + 0.662C Da (F 5j5 = 46.7, p = 0.000338), with 
lowest AIC value of 54.64, was the best one to describe the contri- 
bution of fitness components to the overall fitness cost. In this model, 
Criv Cp ea Cr c , Cri and Qj a are five fitness component costs estimated 
using egg hatching rate, fecundity, emergence rate, larval survival 



rate and developmental duration of adult, and their importance 
decreased according to the significance of their coefficients 
(Table 1). We reported the scatter plots of predicted overall fitness 
cost against observed data in Fig. S3. 

A multiple linear regression model Log 10 -R r = 0.502 + 0.032Cf ec 
+ 0.026Q& - 0.063C Wp + 0.099C R6p (F 4)6 = 10.54, p = 0.007) was 
established to describe the relationship between fitness component 
costs and development of resistance. The significance of coefficients 
and the scatter plot of predicted resistance against observed data are 
reported in Table 2 and Fig. S4, respectively. The variables used to 
predict resistance levels were similar to those predicting overall fit- 
ness cost (described above). Among these four explanatory variables, 
the cost in fecundity Cp eo was the most important, followed by the 
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Figure 4 | Developmental duration (days) of larva (A), female pupa (B), male pupa (C), female adult (D), male adult (E), pupal weight (F) and fecundity 
(G). Bars indicated the standard error (SE). Different letters on the bars indicate significantly differences among the 11 strains (ANOVA followed by 
Tukey's HSD posthoc test, P < 0.05 level). Data and statistics to produce this figure were listed in Table S9. 



cost in copulation rate, C Rc , in pupal weight, C Wp and in survival rate 
from the 6 th instar larvae to pupae, Cr 6p . 

Discussion 

In the present study, the comparison of fitness in 10 resistant H. 
armigera strains with that in a susceptible strain indicated a logistic 
increase in overall fitness costs associated with the log transformed 
resistance level to the Bt toxin CrylAc. Using multiple regressions, 
we established models that predict overall fitness cost and resistance 
level with fitness costs detected in key life history traits. These models 
would enable researchers and extension specialists to estimate overall 
fitness cost or resistance level in the field populations using several 
key life historical traits. Our results will help modelers optimizing 



models that simulate selection of resistance to CrylAc in H. armi- 
gera; this would ultimately help in assessing efficacy of resistance 
management strategies (e.g. for refuge strategies, which are import- 
ant for sustainable use of Bt cotton). 

In the well-known simulation models, on which the current widely 
adopted resistance management refuge strategy is rooted, the fitness 
cost was often treated as a fixed one 10 . Although the fitness cost had 
been found to vary in different environments and higher costs were 
linked with higher resistance levels, the relationship between fitness 
cost and resistance level has not yet been well characterized. In the 
present study, the overall fitness cost seemed to increase in a non- 
linear manner with the Log 10 i?r, which described the relationship 
better than any linear models. The logistic model suggested that 
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Table 1 | Major influential fitness component costs to overall fitness cost and regression coefficients 
Fitness component cost 

(Intercept) 

Crh 2 

Qec 

Cm 

*indicates significant at significant level of 0.05. 
**indicates significant at significant level of 0.01 . 
***indicates significant at significant level of 0.00 1 . 



Coefficient 


Estimated coefficient 


Std. Error 


t value 


Pr(>|t|) 


Significance 


a 


-1.590681 


1.376196 


-1.156 


0.299968 




Pi 


0.030426 


0.003947 


7.708 


0.000587 


*** 


P 2 


0.258139 


0.044570 


5.792 


0.002161 


** 


P 3 


0.617304 


0.1 16802 


5.285 


0.003232 


** 


P4 


0.617594 


0.132943 


4.646 


0.005605 


** 


Ps 


0.662258 


0.185108 


3.578 


0.015914 


* 



the increase in fitness cost of the resistance to CrylAc may asymp- 
totically reach a maximum of 24.47% (±5.22). This is actually quite 
similar to the standard value of 25% suggested for modelers in 2009 
by Gassmann et al. 10 . The logistic model was not only statistically 
significant but also biologically meaningful. 

To explain the logistic increase of the overall fitness cost with the 
Logio-Rr, we hypothesized that development of resistance to CrylAc 
might occur throughout three distinct phases in H. armigera. In the 
first phase, we assumed that only physiological adaptations (a com- 
mon phenomenon in arthropods to any hazard that causes no obvi- 
ous fitness costs 13 ) were involved and the fitness cost increased very 
slowly. In a second phase, we assumed that selection of genetic muta- 
tions conferring high level of resistance was involved. The regulation 
of resistance genes, such as over-expression of detoxification 
enzymes 14 21 , change in hormones induced by detoxification 
enzymes, and structural modification of targeted genes 14,17,22,23 , has 
been considered the most important factor producing fitness cost. In 
this phase, fitness cost increased quickly with an increasing number 
of mutated genes because of the pleiotropic effects of these genes. In a 
third phase, the fitness cost reached a plateau. This might result from 
new mutations providing higher resistance level but fewer pleiotro- 
pic effects, or from compensatory mutations (other mutations) that 
might have been selected. Alternatively, new mutations conferring 
higher resistance levels but lower fitness cost could have been 
selected and have taken over resistance alleles that caused high fitness 
costs. 

The logistic increase of overall fitness cost with the Log 10 Rr was 
observed for H. armigera strains selected on artificial diet with added 
toxin in laboratory, rather than on cotton plants in the field, or 
collected directly in the field. Thus the logistic relationship and the 
above explanations might have some limitations and should be cau- 
tious to apply to the field evolved cases. In the laboratory, the muta- 
tions conferring higher resistance were always selected no matter 
how high the accompanied fitness cost was. However, if the muta- 
tions conferring no or lower fitness cost exist, they would have always 
been selected firstly in the field, because the fitness cost may counter- 
balance the benefit conferred by resistant mutation on refuge crops. 
This has been demonstrated in the field-evolved cost free resistant 
populations of the fall armyworm Spodopterafrugiperda 24 , the west- 
ern corn rootworm, Diabrotica virgifera virgifera 25 and the sugarcane 
borer Diatraea saccharalis 26 on Bt maize. The D. virgifera virgifera 
have no fitness cost even for laboratory selected strains 27 . These three 



species evolved resistance in the field might be lack of mutations with 
fitness costs. This reflected some intrinsic diference between D. vir- 
gifera virgifera and H. armigera or difference in complicated inter- 
actions between pest and host exhibited by corn pests and cotton 
pests. For those species with fitness cost being widely identified in 
laboratory, fitness cost must have been playing an importance role in 
delaying resistance, such as in H. armigera. 

In our study, we have considered the overall fitness cost with the 
actual time elapsed since strains were established in the laboratory 
(i.e. number of generations kept on a constant dose of CrylAc toxin 
since resistant strains were established) to detect potential compens- 
atory effects. However, in our study, the elapsed time since establish- 
ment was significantly related to resistance level. Therefore we could 
not rule out the influence of development of resistance in the 
observed fitness cost increase when testing for compensatory muta- 
tion effects. It would be valuable to collect fitness cost data chronolo- 
gically while the resistance strains are kept on a constant dose of 
toxin, with the ideal situation of keeping resistance level constant. 
Such a study would also help clarify whether a cost-free resistance to 
Bt toxin CrylAc may evolve. In terms of practical applications, such 
possibility of cost-free resistance would prompt re-thinking of the 
refuge method to maintain satisfactory management of resistance in 
targeted pests (H. armigera in Bt cotton in the present case). 

The logistic model was established using strains showing mod- 
erate resistance levels (LF and XJ series) and one showing high res- 
istance levels (XX series). Therefore, our findings would need to be 
confirmed through additional studies, notably using strains that 
would show higher resistance levels in the LF and XJ series. In addi- 
tion, using pooled data from three series of resistant strains might not 
have been optimal conditions for assessing how fitness cost can vary 
with resistance level, notably because the three series of strains might 
have showed differences in genetic background (owing to field popu- 
lations collected from three different locations). However, it has been 
demonstrated that there were four distinct H. armigera geotypes 
(geographical genotypes) in China 28,29 , and the three H. armigera 
strain series used for our study all came from the same genotype. 
Therefore, the three strain series did share similar genetic back- 
grounds, and pooling data was valid in our study. Genetic deviations 
could exist among different populations; however, such differences 
could not have been enough to cause the important differences in life 
history traits that we observed among the populations. This conclu- 
sion is also consistent with the lack of significant difference in linear 



Table 2 Major influential fitness components costs to resistance ratio and 
Fitness component cost Coefficient Estimated coefficient 


coefficients 
Std. Error 


t value 


Pr(>|t|) 


Significance 


(Intercept) a 


0.502961 


0.264468 


1.902 


0.10591 




Cr6p Pi 


0.098970 


0.047561 


2.081 


0.08263 




CWp P2 


-0.062893 


0.025234 


-2.492 


0.04701 


* 


c Fec p 3 


0.03171 1 


0.007606 


4.169 


0.00589 


** 


Crc P4 


0.026279 


0.008480 


3.099 


0.021 15 
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regression slopes between LF and XJ series. Therefore, the costs 
detected in the various strains can be attributed confidently to the 
development of resistance to CrylAc. 

A meta-analysis suggested that fitness costs were more frequently 
detected in selection experiments than in fitness component com- 
parisons 10 . The present study provided an explanation by analyzing 
the correlation between overall fitness cost and each fitness compon- 
ent cost. With multiple regression analyses, we demonstrated that the 
egg hatching rate (squared), fecundity, emergence rate, larval sur- 
vival rate, and developmental duration of adults (ranked in decreas- 
ing importance) affected the overall fitness cost. Although costs in 
developmental duration of larvae and pupae were highly correlated 
with the overall fitness cost (pairwise correlation coefficients, Fig. 
S2CL and EL), both were discarded through multiple regressions. 
This might be because (1) the developmental duration of larvae 
varied among strains in LF series, whereas it did not in XJ series 
(Table S9) and (2) the developmental duration of pupae did not 
change much in strains showing low resistance levels, but it did 
increase in strains showing the highest resistance levels in each series 
(Table S9). The strains with the highest resistance level in XJ and LF 
series also had the shortest elapsed time for establishment (Fig. 1). 
This implied that the cost in developmental duration of pupae might 
have been compensated in those strains with lower resistance levels; 
therefore, this trait should be considered in future studies aiming to 
understand the compensation of fitness costs with development of 
resistance to toxins in arthropods, or at least in the case of H. 
armigera. 

Overall, a low probability of detecting fitness cost when comparing 
life history traits may be due to (i) an absence of cost related to 
CrylAc resistance in various life history traits in H. armigera (e.g. 
sex ratio), (ii) a quick compensation of the cost in some traits (e.g. 
developmental duration of pupa), and to (iii) random variations in 
traits evaluated. 

We identified the life history traits linked to the development of 
resistance to CrylAc in H. armigera. Both fecundity and larval sur- 
vival were identified as important traits influencing overall fitness 
cost related to the resistance to CrylAc. The survival from 6 th instar 
larva to pupal stage was actually more sensitive than survival of the 
whole larval stage when assessing development of resistance to the Bt 
toxin. Pupa is the life stage undergoing important physiological and 
morphological transformation; it can therefore be the crucial stage to 
reveal disadvantages linked to the selection of mutations conferring 
resistance to Bt toxin. Because CrylAc toxin primarily targets the 
midgut of the larval stage of Lepidoptera, and reduction or loss of 
binding of Bt toxin to midgut membrane receptors is the most com- 
mon general mechanism of resistance 30 " 32 , it has been hypothesized 
that mutations conferring resistance to Bt toxins might cause fitness 
costs by reduction of nutrient uptake in guts of the Lepidopteran 
pests (e.g. by increasing permeability of the gut membrane to toxic 
phytochemicals or pathogens, or by impairing protein digestion) 10 . 
Our findings suggest that such mutations may be particularly det- 
rimental when the larvae pupate. 

Methods 

Insect strains. Eleven H. armigera strains were reared under laboratory conditions 
(27 ± 2°C, 75 ± 10% RH, 14L:10D). According to the origin of the strains, they were 
grouped into 3 series: (i) the XX series included the 96S and BtR strains (ii) the LF 
series included the LF5, LF10, LF20, LF30 and LF60 strains, and (iii) the XJ series 
included the XJF, XJ1, XJ5 and XJ10 strains (Fig. 1). The original populations of the 
three series of strains were all collected from Yellow River Region and were subject to 
various screening programs to produce strains with different resistance levels (Fig. 1 
and supplementary text). All bioassays described below were carried out in similar 
conditions. 

Bioassay of resistance to CrylAc toxin. The H. armigera strains were reared and 
assessed simultaneously for susceptibility to CrylAc using the methods reported by 
Liang et al. 12 ' 33 . In brief, the susceptibility was evaluated by feeding H. armigera larvae 
with artificial diets containing various concentrations of CrylAc toxin. Four-day-old 
larvae were placed individually on Cry 1 Ac-contaminated artificial diet in plastic wells 



(depth: 1,5 cm, vol. 3 ml) in 24- well plates. The plates were covered with a plastic lid 
to prevent escape and 96 larvae were tested per CrylAc concentration. After seven 
days, we recorded the mortality rate and measured body mass of individuals still alive. 
Both dead larvae and those with a body mass of less than 5 mg were recorded as 
dead 34 . For each strain, the Lethal Concentration 50 (LC 50 ) value was determined 
using a probit analysis 35 and the resistance level to CrylAc indicated by the resistance 
ratio was calculated with LC 50 of the tested strains divided by the LC 50 of the 
susceptible strain. 

Assessment of fitness cost. Thirteen fitness components {life history parameters) 
were obtained with life table techniques {Supplementar text). The components 
considered were (i) survival rate from the 1 st to 6 th instar larvae, (ii) survival rate from 
the 6 th instar larvae to pupae, (iii) duration of larval development, (iv) pupal weight, 
(v) developmental duration of female pupa, (vi) developmental duration of male 
pupa, (vii) emergence rate, (viii) proportion of females (or sex ratio), (ix) copulation 
rate, (x) developmental duration of female adult and (xi) of male adult, (xii) fecundity, 
and (xiii) hatching rate. 

Each fitness component cost, Q, for each strain was calculated as: 

Ci = [(F R -Fs)/Fs] x 100% for developmental durations 



or 

Q — [(Fs~Fr)/Fs] x 100% for all other fitness components 

where F s and F R are the mean values of fitness components for susceptible and 
resistant strains, respectively 10 . 

For each strain, the intrinsic rate of population increase, r m , was calculated basing 
on the life table data with the methods described by Bird and Akhurst 36 : 

r m = (]nRo)/T 

where T is the developmental time from egg to adult eclosion. 

Ro = (n xl e xI a )/2 

where n is the mean number of eggs per female, I e is the proportion of fertile eggs, l a 
is the proportion of reproductive adults produced from the 1 st instar larvae, and 2 is 
the sex ratio coefficient for H. armigera. 

The overall fitness cost C was then calculated as: 

C=[(r mS -r mR )/r mS ]xl00%, 

where r mS and r mR are the intrinsic rate of population increase for susceptible and 
resistant strains, respectively. 

The Q and C were calculated with reference to either the susceptible strain 96S or 
the relative susceptible strains in each series (i.e., LF5 and XJF for LF and XJ series, 
respectively), and calculated negative values (if any) were converted to 0. 

Statistical analysis. The survival rates were calculated using the Bayesian method and 
compared among strains with Monte Carlo approximations 37 (Supplementary text). 
The means of developmental duration, weight and fecundity were compared among 
strains using an ANOVA. The resistance ratio was loglO transformed to normalize 
the data distribution before regression analysis. 

To assess how the fitness cost varied with the development of resistance, a series of 
linear and nonlinear regressions were performed for each series of strains separately 
and for pooled data (Supplementary text). Multiple linear regressions were used to 
analyze how each fitness component cost influenced the overall fitness cost and 
resistance level (Supplementary text). 

The above analysis was carried out for overall fitness costs calculated with reference 
to either the susceptible strain 96S or the relative susceptible strains in each series (i.e., 
LF5 and XJF for LF and XJ series, respectively). Only the former produced satisfied 
results (i.e., statistically significant results) and reported here. All statistics were 
performed with R 38 . 
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